THE IATERAL/DIKECTIONAL STABILITY CHARACTERISTICS 
OF A FOUR-PROPELLER TILT -WING V/STOL MODEL IN 
LOW-SPEED STEEP DESCENT 
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SUMMARY 

Lateral -directional dynamic stability derivatives are presented 
for a 0.1-scale model of the XC-142A tilt-wing transport. Thj stability 
derivatives were determined from experiments conducted in tLo Princeton 
University Dynamic Track with a dynamically similar model. The tests 
involved various descending flight conditions achieved at constant 
speed and wing incidence by varying the vehicle angle of attack. The 
propeller blade angle and the speed were also changed in the steepest 
descent case. 

The experimental data were analyzed assuming that the dynamic 
motions of the vehicle may be described by linearized equations, with 
the lateral-directional characteristics of the full-scale aircraft 
also presented and discussed. Results from this experimental investi- 
gation indicated that the full-scale aircraft would have a stable 
lateral -directional motion in level flight, with the dynamic motion 
becoming less stable as the descent angle was increased. No improve- 
ment. in the dynamics was noted when the propeller blade angle was 
reduced daring the steepest descent, although a subsequent increase 
in airspeed caused a further degradation in the lateral-directional 
stabi] tty which was characterized by an unstable Dutch -roll oscillation. 
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INTRODUCTION 


During the past decade, there has ^een an ever-increasing interest 
in the V/STOL type of aircraft, not only because of their unique and 
intriguing addition to the field of aerodynamics, but also because of 
the potential offered by these aircraft in a commercial and military 
role. An assortment of V/STOL configurations have been tunnel tested, 
including tilt-prop, compound, jet-lift, tilt-wing, and many other 
concepts (refs. 1 through 9) • SUch configurations offering promise, 
either have been built and flight tested, or are currently involved in 
a flight research program (refs. 10 through 18 ). Most of this work, 
though, has been of a gen ral exploratory nature, directed towards 
either pure aerodynamics or establishing the overall feasibility of the 
various concepts. Investigations to document the classical stability 
derivatives, which are a prime ingredient affecting the so-called 
handling qualities of the vehicle, have been limited. 

To make this stability derivative data available, a number of 
V/STOL configurations have been tunnel tested in the Pr'nceton University 
Dynamic Model Track (refs. 5 to 7, and 19 to 21). One of the more recent 
vehicles tested has been a 0.1-scale model of the four -propeller XC-1U2A 
tilt-wing transport. Experimental results have included the first quan- 
titative information published on the lateral-directional dynamic 
derivatives of a tilt-wing V/STOL aircraft at low forward speed (see 
ref. 19 ) and the longitudinal dynamic stability characteristics (see 
ref. 20). 

The purpose of this study was to extend the information presented 



in the previous Princeton University reports dealing with the XC-lk2A. 
Direct quantitative measurements were made of the lateral -directional 
tra n sient response time histories for various trim conditions simulating 
low-speed, descending flight. These data were taken at five test con- 
ditions representing various combinations of pitch attitude (0)> for— 
velocity (U f ), and propeller blade angle (p). The wing incidence and 
wing flap angle were fixed at 40° and 60°, respectively, throughout the 



entire investigation. The resul’ts presented in this report thus 
represent the first documentation of the stability derivatives for a 
tilt-wing V/STOL transport in descending flight. 

The servo -analysis techniques used to determine the stability 
derivatives from the data available are described in appendix 1. 
Conversion of the data to full-scale values is discussed in appendix II* 
with the results shown in figure 1. Additional information regarding 
the static characteristics and control effectiveness for the configu- 
ration tested, are presented in reference 22. 



3 


Mi' 




DESCRIPTION OF APPARATUS AND EXPERIMENTS 

TEST FACILITY 

Tie Princeton University Dynamic Model Track Facility was designed 
explicitly for the study of the dynamic stability and control of heli- 
copter and V/STOL models for speeds ranging from hover through transi- 
tion. Integral components of the test facility include: a 7?0-foot 

track, servo-drive carriage, model mounts, measuring transducers and 
recording equipment, all of which are located in a building with a test 
cross section of 50 by 50 feet. The dynamic carriage, which can follow 
the longitudinal velocity excursions of the model, has an acceleration 
potential of 0.6g and a maximum speed of 40 feet per second. A detailed 
description of the facility and testing techniques employed may be 
found in reference 25 . 

Two of the various methods used to mount models to the carriage 
permit separate measurement of the classical longitudinal and lateral- 
directional degrees of dynamic motion. The longitudinal mount, shown 
in figure 2, permits horizontal and vertical motions of the model 
relative to the carriage and allows the model to rotate in the plane 
determined by these directions. The horizontal motion is sensed and 
used to command the carriage to follow the model in a closed-loop 
fashion. Similarly, the vertical displacement of the model comnands 
the boom to folio# in the vertical direction. The lateral-directional 
mount which was used for this study is shown in figure 3* This mount 
permits relative motion between the model support linkage and the 
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lateral servo-driven carriage to be sensed and used to position the 
lateral carriage along the lateral boom. A maximum sideward excursion 
of 8 feet is permitted by this arrangement. The yaw degree -of -freedom 
is provided by a pivot mounting which allows angular rotation of the 
vertical tube supporting the model relative to the lateral servo- 
driven carriage. Roll freedom is achieved by a pivot mounting located 
within the fuselage of the model itself which permits the angular 
motion in roll relative to the fixed vertical support tube. A schematic 
drawing of the lateral -directional mount is shown in figure 4. Overall, 
the model support and gimbal system allows particular degrees of angular 
freedom to be selected. Those not under investigation are locked by a 
braking system which also serves to arrest the model motions at the end 
of a test run. The gimbal and support systems also serve as references 
for measurement of the model motions. This system is similar to the 
one vised in reference 19, with the roll and yaw axes fixed to the 
lateral error linkage; the roll axis yaws but does not pitch with the 
model, while the pitch axis remains fixed to the model. For the pre- 
sent study, the descent conditions were simulated by altering the pitch 
attitude of the model with respect to the carriage. This adjustment 
results in an increased aircraft angle of attack since the carriage 
velocity is always horizontal. The exact expressions for the vari- 
ables and the appropriate equations of motion for this axis system 
are presented in appendix I. 

Hie dynamic experiments conducted during this study were, for 
the most part, two degree-of-freedora motions achieved by use of the 
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lateral -directional mount. These tests included single degree-of- 
freedom in yaw (t) and two degree -of -freedom in roll and yaw (0 - ty) 
angular motions, and roll angle and sideward velocity (0 - v^) motions 
for selected combinations of pitch attitude .(9), forward velocity (U^.) 
and propeller blade angle (p). Three degree -of -freedom cases 
(0-4 - Vf) were also conducted to analyze the complete lateral- 
directional motions of the aircraft. 

The Princeton Model Track facility can also be used to measure 
static stability derivatives. The model i r mounted rigidly to the 
carriage permitting the forces and moments acting on the model to be 
measured with strain gauges. Continuous changes in the variables of 
interest can be achieved by programming the model or carriage motions. 
Although si m i l a r to wind-tunnel testing, a uniform air velocity, free 
from turbulence, is offered in the JO by JO foot test section. Precise 
control of speed is available over a wide range, which includes not 
only forward flight but hover and rearward flight as well. 

MODEL 

The 0.10 scale dynamically similar model constructed for this 
experiment is shown in figure 5; drawings of the general arrangement 
and of the airfoil section of the model are shown in figures 6 and 7* 
The model is the one used for the lateral-directional study of 
reference 19, with several modifications. These modifications of the 
model, and descriptions of the lateral-directional control system, 
geometric characteristics of the model propellers, ailerons, and other 
pertinent information are presented in reference 22. 




INSTRUMENTATION 


The basic test instrumentation, employing telemetering and magnetic 
tape recordings, was similar to that used in reference 19* However, 
several additions and refinements were made -to the model instrumenta- 
tion: these changes are also noted in the data report relating to 

this experiment (ref. 22). 





EXPERIMENTAL RESULTS 


STATIC DATA 

Measurements of the horizontal force (parallel to the free -stream) 
and the vertical force (perpendicular to the free-stream) were made to 
determine the descent conditions realized during this program. All 
data were taken in the forward flight mode for a constant wing inci- 
dence, i = 40°, and constant flap deflection, 8 f = 60°. 

Generally, the aerodynamic forces are functions of the forward 
velocity (U^), vertical velocity (W^), pitch angle (0), blade collec- 
tive pitch angle (p), and the control inputs. However, for the XC-142A 
model, there was no cyclic pitch control and, the ailerons, differential 
collective, etc., were assumed to affect only the lateral-directional 
equilibrium. Therefore, for a space-fixed axis system, the force 
equations for steady level trimmed flight are: 


x f (u f ,w f ,e,p) . o 

Z f( U f' W i* 0 ,P)“ “8 


( 1 ) 


Since a space -fixed axis system is being considered, the descent 
angle (7) is defined by: 


V 

tan 7 - -p (2) 

However, because of the mechanization of the testing apparatus. It 
was not possible to program descent velocities Into the model, that 
Is, V f - 0. 
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Therefore, the descent conditions were simulated hy rotating the model 
about its pitch axis, thus stipulating that the test or true model 
pitch attitude ( 0 ) was precisely the model angle of attack (a), as 
suuwu on figure 8. For each of the test cases studies, this pitch 
attitude remained constant. Force balance considerations were then 
used to define the actual descent angle, by, 

Zf sin 7 = cos 7 (3) 


or 

Xf 

tan 7 = — 

Z f 

Since the horizontal velocity of the model is determined by programming 
the carriage, it is not necessary in the model tests that the horizontal 
forces be in balance. 

It should be noted that the model was constrained to fly parallel 
to the horizon along the space X-axis, (i.e., the test track), thus, 
simulating an inclined horizon as shown on figure 8, rather than an 
actual descending flight path. Since the horizon is so inclined, an 
aerodynamic pitch attitude (© a ) can be defined as the angle between 
this artificial horizon and the body principal X-axis, again noted 
on figure 8* This will be the pitch angle of the aircraft in descending 
flight. 

She space axis longitudinal force (Xf ) was measured directly by 
a strain gauge, thus known for the various combinations of U f , 6, 
and Direct measurement of the apace-exia vertical force (Zf) was 
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not available. To determine Z f , the model was rolled through 0 
degrees as shown in figure 8, and the horizon referenced side force (Y) 
measured. The relationship for obtaining the space vertical force 
is then 


- Y _ Y 
Y sin 0 0 


(for small angle 
approximation ) 


When subftituted into the descent angle equation, this yields: 


7 s tan "y" (U^, 6,3)0 


Samples of the static data used to determine the horizontal and ver- 
tical forces, and thus the descent angles, for each test case are 
shown in figures 9 through 11. A summary of the test conditions. 
Including the descent angles computed by the above equation, and the 
aerodynamic pitch angle (@ a ) are presented in table I, with the model 
geometric and inertia characteristics outlined in table II. The 
details of other pertinent procedures and testing techniques appli- 
cable to this study, regarding the stability augmentation used on the 
model, model trim, and aileron effectiveness, are discussed in 
reference 22. 


DYNAMIC DATA 


Single -degree -of -freedom dynamic responses mere evaluated for 
each of the five cases to determine the yaw damping (l^) and the 
moment of inertia about the yarn axis (I ). For these runs, the model 

m 








10 


was locked in roll. Mechanical springs were attached between the model 
and support tube to provide a restoring moment about the yaw axis. The 
equation of motion for the rigid-body oscillation of the model in yaw, 
’.1th mechanical springs, aerodynamic damping, and mounting friction is: 


ij + ^ + frictionjj + ^ t - 0 


For these experiments, the model mounting friction was determined from 

Sn 

power-off tests . The mechanical' spring contribution, ~ = K^, was 
calibrated separately and found to be (-31*7 ft-lb/rad) . Spring- 
restrained single -degree -of -freedom runs were also performed with the 
model motors off and the carriage velocity equal to zero (U^ = rpm= 0) 
to determine the model yaw moment of Inertia. The yaw damping runs 
were conducted by releasing the model, with propellers running, from 
an initial yaw angle offset. An example of the resulting time history 
plus one for the model motors off are shown in figure 12. These time 
histories, along with the knowledge of the spring constant and moment 
of Inertia, provided the total aerodynamic plus mechanical spring 
damping. The mechanical damping, as determined from model motors 
off data, vas simply subtracted out, resulting In the true aerodynamic 
damping for the first case. Since the model exhibited stable direc- 
tional stability In the first test condition, the mechanical springs 
were not used during the remaining descent cases. It should be noted 
that no single -degree -of -freedom runs mere made to determine the roll 


damping and roll mo m ent of Inertia. As an approximation, therefore, 
the ratio of the moment of Inertia In yaw far this experiment to that 


11 


of reference 19 was used to determine the roll moment of inertia. 

Furthermore, the roll damping moment for the level flight case of 
this investigation was assumed to be equal to the value determined 
in the level case (i w = 30 °) of reference 19 , resulting in the slightly 
lower ratio of roll damping to inertia (hereafter simply the roll 

angular rate damping). A summary of the single -degree -of -freedom data ; 

for the test runs analyzed for each case is presented in table III. * 

j 

Multiple Degree of Freedom 

For each of the descent cases, lateral -directional transient 
response measurements were conducted for the following degrees of 
freedom: two degrees of freedom in roll and yaw (0 - ♦); two degrees 

of freedom in roll and sideward velocity (0 - v^); three degrees of 
freedom in roll, yaw and sideward velocity (0 - ♦ - v^). The equations 
of motion for each of these conditions sure presented in appendix I. 

No mechanical springs were used for any of the mult iple -degree -of - 
freedom tests, and for the data presented, only pitc hing motions 
were stability augmented. 

The procedure followed during one of these mult iple -degree -of - 
freedom tests was to bring the model motors up to speed with the 
model locked, relative to the boom, on the stationary carriage. The 
carriage was accelerated to the trim speed, and timers were used to 
unlock the mechanical restraints of the desired degrees of freedom. 

The model was then permitted to fly freely in any desired combination 
of degrees of freedom* For some of the runs, either a roll or ding 
offset, or control deflections were used as an input to excite the 
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model's response, though in some cases a small random disturbance was 
sufficient to start the motion due to the unstable nature of the 
dynamics. 

Sample results of the lateral-directional dynamic tests for the 
five cases, which include level flight, are presented in figures 12 
through 31 as time histories of the transient response of the model. 

The period of the motion and the damping ratio for each of the multiple- 
degree -of -freedom runs analyzed is presented for each case in table IV, 
with the corresponding roots of the model transient responses sumnarized 
in figure 32. 
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FUIiL-SCALE AIRCRAFT DYNAMICS 

Hie analysis of the experimental results, as discussed in appendix I, 
along with the application of the scale factors discussed in appendix II, 
permitted the full-scale, body axes stability derivatives to be deter- 
mined. These stability derivatives are presented in figure 1 as u 
function of descent angle, •j'he dashed curves indicate the probable 
trends of the derivatives, as determined from test cues 1,2, and 3» 

The effects of a change in the propeller blade angle (0), case 4, and 
subsequent change in the forward velocity (U^.), case e >. are also shown 
in figure 1, by the rectangular and triangular symbols respectively. 

Also noted in figure 1 are the results from reference 23, which will be 
discussed later as applied to the level flight case. The data resulting 
from descent angle changes, cases 1,2, and 3, will be discussed separately 
from the data representing aerodynamic changes, cases 4 and 3* 

DESCENDING FLIGHT RESULTS 

The descending flight results were those obtained from cases 1 
through 3. TVo descent angle obtained by varying fusalage angle of 
attack was the only variable for these three cases and were: 

°° (level), -5°, -11°, respectively. From the data analysed, the 
body axis stability derivatives were computed and the following nay 
be noted for each. 

Roll Damping (Lp) 

A ma the m at ical approximation was rsployed to determine the rail 
angular rate damping of the model, mm a taction of descent angle, 
because no slngle-dagree-of -freedom roll date was available. It was 
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not known whether the roll damping would increase or decrease as the 
descent angle steepened, therefore, variations in both directions were 
analyzed. The results showed that an increase in the space axis roll 
damping (as would have been measured from the track data) caused the 
body axis values to increase unrealistically for the higher descent 
angles. On the other hand, a reduction in the space axis roll damping, 
for higher descent angles, resulted in a ~..igbt Increase in the final 
body axis value, as shown in figure 1, cases 1 to 3* This trend in the 
roll rate damping was acceptable because these values, along with the 
other stability derivatives involved, were consistent with the dynamic 
motions measured for these descent cases. Such a variation in the roll 
damping was also considered reasonable because it can be attributed to 
an Increase in lift on the wlng3, at the higher angles of attack associ- 
ated with these steeper descent angles, a condition that would be 
expected prior to vlng stall. 

Yaw Moment !Xie to Boll Rate (Up) 

This cross derivative is negative - the sign generally associated 
with conventional aircraft (ref. 2t) - resulting from the unsyametric 
lift distribution on the wings while the aircraft is rolling. The 
magnitude remains som e w h at constant, though falls off slightly at the 
higher descent angle. This reduction may be attributed to a smaller 
tall contribution and/or unsymmetrical stall effects on a tilt-wing 
configuration that Is rolling during descending flight. 
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Yaw Damping (N r ) 

This stability derivative is stable for level flight and becomes 
less stable as the descent angle increases. For the largest descent 
angle, in fact, th« damping becomes very slightly unstable (positive) 
vMch arises, most likely, from interference on the vertical tail 
particularly at the low speeds being investigated. 

Roll Moment Due to Yaw Rate (L,.) 

The rolling moment due to yaw rate is another cross derivative 
which has a sign (positive) that is typical for conventional aircraft 
(ref. 24). The value of the derivative decreases for the first descent 
case (-5°) but increases positively for the higher descent case. The 
explanation of this trend appears somewhat difficult and is probably 
the result of the complex airflow interplay between the wing and tall 
surfaces . 

Dihedral Effect (L v ) 

This stability derivative exhibited a stable (negative) tendency 
for all the descent cases analyzed. Its value is large for the level 
flight case - comparable to that found in reference 19 and becomes less 
stable as the descent angle increases This trend may be caused by an 
onset of stall on the advancing wing, permitting the opposite wing's 
lift to become more effective with descent angle. 

Directional Stability (N y ) 

The sign of this stability derivative Indicates the aircraft is 
direct ionally stable (positive) over the range of test conditions* 

The magnitude appears to remain constant with a slight increase in 
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stability noted for the larger descent angles. The value of N y for 
the level flight case is in good agreement with the value calculated 
from the approximation based on the fact that the free -stream velocity 
on the vertical tail is the primary factor contributing to this deriva- 
tive. From reference 24, assuming negligible sidewash, and a vertical 
tail efficiency factor of one, the directional stability can be written 
as: 


N v = 


qSfpArpliji 

“i^oT 


and was computed using the following 

A^, = 2 per radian 

I z = 270,000 slug-feet squared 

Z T » 23.5 feet 

Sip = 130 square feet 

U 0 =57.8 feet per second 

resulting in 

N y » 0.0016 per foot-second 
Sideforce Due to Lateral Velocity (Y v ) 

This stability derivative was found to be small and alsc remained 
fairly constant throughout the range of test conditions. 

For the above descent cases, the data indicated that the dynamic 
motion of the full-scale aircraft in level flight was stable. The 
rolling mode would be convergent, with a time to one-half amplitude 
of 0.91 second , while the spiral mode would also be stable having 
a time to one-half amplitude of 49*5 seconds. The Bitch-roll 
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oscillation would be approximately neutral with a period of 8.4 seconds 
and a time to one-half amplitude of I 38 .O seconds. The results for the 
second and third descent cases indicate the dynamic motion of the full- 
scale aircraft was similar. That is, both cases exhibited a spiral 
divergence, a rolling convergence, and a stable Dutch-roll oscillation. 
Specifically, for the 5° descent (case 2), the period of the Dutch-roll 
mode would be 7*9 seconds and a time to one-half amplitude of 53*0 
seconds. The rolling convergence would have a time to one-half ampli- 
tude of 0.9 second, while the time to double amplitude of the divergent 
spiral mode would be 8.4 seconds. For the 11° descent (case 3) the data 
indicated the Dutch-roll period would be 8.4 seconds with a time to one- 
half amplitude of 8.1 seconds. The rolling convergence would have a 
time to one -half amplitude of 0.9 second, while the time to double 
amplitude for the spiral divergence would be 2.4 seconds. 

AERODYNAMIC CHANGES 

The test conditions of eases 4 and 3 were added to determine the 
effects of very limited aerodynamic changes for a particular descent 
case. Specif ically, a reduction in the collective blade angle (0), 
from 13 . 5 ° to 11 . 3 °, was he only change made for case 4; further, at 
this new blade angle setting, the model test velocity (U f ) for case 5 
was increased from approximately 18.5 feet per second to 24 feet per 
second. The equivalent full-scale velocity for these two values are 
58.5 fps (34.5 knots) and 76 fps (45 knots), respectively. 

The stability derivatives calculated for the last two cases are 
compared to those of the descent cases on figure 1. Generally, the 
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results indicate that very little change, if any, was realized for the 
reduced blade angle setting case. Increasing the velocity shows that 
the effect on the stability derivatives was in the same direction as 
compared to case 4, though, for the most part, resulting in a much 
larger magnitude change. The only exceptions were those trends 
indicated for tte dihedral effect and the sideforce derivative which 
were opposite to those realized in case 4. 

The forward velocity effects, particularly regarding the side- 
velocity derivatives, are at best very difficult to explain. This is 
so because of the unknown interaction between the body, wing, and tail 
airloads, as affected by the free-stream and/or slipstream, particularly 
for a descending tilt-wing configuration, itor example, a variation in 
lift, which may explain one stability derivative change would not 
necessarily explain a change in a different stability derivative. 
Generally, it can be said that the decrease in N y Is caused by a 
wake turbulence or interruption of the flow at the vertical tail - 
for the steep descent cases. An increase in the sideforce derivative 
can be attributed to fuselage impingement, while the larger negative 
value of L y results from the change in lift experienced on each 
wing. 

The dynamic motion of the full-scale aircraft computed for case 4 
was almost identical to that of case 3, Apparently, as discussed above, 
a change in the blade collective pitch setting alone did not greatly 
alter the dynamics. The model results indicate that the full-scale 
aircraft would exhibit a stable Dutch-roll oscillation with a period 
of 8.1 seconds and a time to on e - h al f amplitude of 9.1 seconds. The 


19 


rolling convergence would, have a time to one-half amplitude of 1.0 
second, while the spiral mode would "be divergent, with a time to double 
amplitude of 2.6 seconds. The model results for case 5 indicate that 
the full-scale aircraft v r \ld exhibit quite a different dynamic motion. 
Specifically, an unstable Dutch-roll type oscillation was •"'■’••played, 
one that would have a period of 8.2 seconds and a time to double 
amplitude of 49.0 seconds. The rolling mode remained convergent, and 
would have a time to one -half aufplitude of 0.9 second, while the 
spiral mode would be less divergent with a time to double amplitude 
of 9*6 seconds. 

The results for the level flight case are, for the most part, in 
good agreement with the level flight results of reference 19. The only 
stability derivative that differs noticeably is the yaw moment due to 
roll rate (Np). The value found in the present study is of much larger 
magnitude and different in sign, that normally exhibited by conven- 
tional aircraft (i.e., negative). It is this term that apparently 
contributes to the slightly stable spiral mode, noted herein, another 
result that differed with the level flight case of reference 19. 

These differences would be attributed to: 

(a) The slight alterations in the testing procedure used during 
these investigations are noted in reference 22. 

(b) The configuration of the model used in the present study 
to achieve level flight, particularly iy * 40°, as compared to the 
iy » 30° configuration vised in reference 19. 
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The overall descent results appear to lend themselves to the dis- 
cussion on descending flight in transition found in reference 2. 

Similar to the descent results noted in this reference, no wing 
buffet was encountered throughout the descent angle range of the 
present study and, by this fact, the feasibility of much higher des- 
cent angles is not known. It is believed, however, that any additional 
descent angle capability would be slight, because of the proximity of 
wing stall and the onset of poorer dynamics indicated by case 5* 

It is also interesting to note that for direct data comparisons 
the computed stability derivati „ es presented herein compare favorably 
to those reported in reference 25. For the appropriate wing angle case 
(i w = 40°) that is comparable to the level case of the present result, 
all the stab? licy derivatives are in close agreement with the exception 
of N v and L^. These differences are evident in figure 1 and may be 
attributed to the variation in the testing techniques between the two 
studies. In addition, the trends in the roll and yaw damping stability 
derivatives as a function of descent angle, noted in reference 25 , are 
in good agreement with the present results, again when the testing 
procedures are considered. Specifically, the descent conditions of 
reference 25 were achieved by changing the propeller rpm for a con- 
stant model angle of attack and wing incidence. This was not the 
procedure used to simulate descent for this study; therefore, it is 
fair to say that the descent cases per se are not directly comparable. 
The test conditions of the present Investigation that are more in 
line with those of reference 25 , are simply those for which the 
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angle of attack was held constant while some other parameter related 
to a change in thrust was varied. Clearly, this type of change is 
represented by cases 3 to 5- Comparison of the stability derivatives 
for these cases with those of reference 25, -that is, yaw damping, roll 
damping, and directional stability, indicate that the results are in 
close agreement both in magnitude and trend. The only exception is 
that the yaw damping results do not agree. Again, this difference 
may be accounted for by the testing technique used, in that the damping 
term found during the forced oscillation tests of reference 25 was 
(N r - Np), from which the yaw damping (N r ) can not be separated. (Note, 
N. results from the yaw moment due to the rate of change of sideslip 
angle . ) 
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CONCLUSIONS 

A O.l-scale model of the XC-142A V/STOL transport was tested in 
the Princeton University Dynamic Model Track while in low-speed level 
and descending flight. The test cases included various descent angles 
at constant speed, and two conditions incorporating aerodynamic changes 
at the steepest descent. On the basis of the model results, the full- 
scale aircraft in level flight would be dynamically stable. A degrada- 
tion in the aircraft's dynamics would be experienced as the descent 
angle Is increased, to a point that for an 11° descent, the motion 
would be made up of a rolling convergence, and a lightly damped Dutch- 
roll oscillation and a spiral divergence with a time to double 
amplitude of 2.4 seconds. 

The first aerodynamic variable change, consisting of a propeller 
blade angle change for an 11° descent, did not significantly alter the 
aircraft's dynamic motion, as compared to the original 11° descent 
case. However, for the second aerodynamic variable change, which 
consisted of increasing the airspeed in the reduced blade angle con- 
dition, the lateral-directional dynamics become more unstable, charac- 
terized by an unstable Dutch-roll oscillation with a time to double 
amplitude of 49.0 seconds. 

Generally, the overall results appear consistent with the theory 
used for the analysis and computation of the stability derivatives. 

In addition, the results and trends are in good agreement with the 
limited data available from similar studies previously conducted. 
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APPENDIX I 
ANALYSIS OF DATA 

AXIS SYSTEM AND EQUATIONS OF MOTION 

Model motions were measured with respect to a moving carriage; 
therefore, it is more convenient to use a spue v, -fixed axis system to 
analyze these motions rather than the more conventional stability 
axis system. The X^-axis is located along the direction of the motion 
of the carriage, with the Z^-axis perpendicular to it, positive down- 
ward; the origin of the axis system is placed at the pivot point where 
the model is attached to the vertical link. The general body axis and 
space axis relationship is shown in figure 33 . For the analysis pre- 
sented herein, tj, the angle between the model principal axis and fuse- 
lage reference, is assumed to be zero. Thus, the angle £, shown in 
figure 33 is exactly equal to 0. The order of rotation U6ed to 
relate the body and space axes are as follows: 


( 1 ) 

* 

rotation about 
roll axis 

z f 

(yaw gimbal) to aline Xf 

with 

(2) 

0 

rotation about 
pitch gimbal 

X’ 

(roll gimbal) to aline Y' 

with 

(3) 

0 

rotation about 
aircraft axis 

T 

(pitch gimbal) to aline X* 

’ with 


The resultant aircraft at- rates expressed in terns of the space 
axis rates are: 

p ■ i cos 0 - ♦ cos 0 sin 0 
r - v cos 0 cos 0 + 0 sin 0 
q - 0 + ♦ sin 0 




( 7 ) 
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When 0 Is assumed small, the above relationships reduce the 
equations (2) of reference 22 

• • 

p ** 0 cos 0 - ♦ sin 0 
r = ♦ cos 0 + 0 sin 0 

q a 0 

Differentiating, the angular acceleration expressions are: 

p » 0 cos 0 - ♦ sin 0 


r » f cos 0 + 0 sin 0 


q - 0 


Using the vector equivalence for rates and moments, that Is, 


P - I® 
rWK B 


e-ifc 


the rate expressions become: 

Lg » Iq cos 0 - % sin 0 

%-^eos •4>X 0 slne (11) 

Since the current study deals only with the lateral-dlreotioael stability 
derivatives, the analysis Kill be continued with tbs applicable rolling 


and yewiap 


expressions. Using the first two 






of equations (ll), to solve for the gimbal moments as a function of 
the tody axis moments, the results are: 


I^j = Lg cos 9 + Ng sin 9 


% = % cos 0 - I'B sin 0 


( 12 ) 


The general principal axis equations of motion (refs. 24 and 26) are: 


I x p + qr(I z -I y ) - (r + pq)l xz * Lg 
I z r + pq(l y -I x ) - (p - qr)!^ = % 

Dropping the terms having products of rates and substituting these 
into equations (12) 


(13) 


% “ co3 9 - I^r cos 0 + I z r sin 0 - I^p sin 0 

% » I z r cos 0 - I^p cos 0 - ^p sin 0 + I xz r sin 0 

Finally, substituting for the body-axis angular accelerations and 
combining coefficients of the space-axis angular accelerations: 

cos 2 ® - 21 x 2 sin 0 cos 0 + sin 2 00 

+ [(I 8 - I x )sin 0 cos 0 + (sin 2 ® - cos 2 ®)IxJt 

Iq - [^sin 2 ® ♦ 2Ixm sin ® cos 0 ♦ cos 2 0t 

♦ [(I* - I*) sin O cos 0 ♦ (sin 2 ® - oos 2 ®)lj0 


<*) 


( 15 ) 
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letting, 

I v cos 2 e - 2I V _ Jin 9 cos 0 + I„ sin 2 0 a I' 

X XZi Lt X 

p p 

I x sin 0 + 2I XZ sin 0 cos 0 + I z - cos 0 a I z * (16) 

(l z - I x )8in 0 cos 0 + (sin 2 0 - cos 2 0)I xz a 


the general equations 1'or the roll and yaw gimbal axes are then 


v* ♦ v* • h 

V* + f xzi * »0 


( 17 ) 


The appropriate gimbal axis moments and cross product of inertia were 
computed for each descent case as a function of the pitch angle 9, 
and are shown below. (Body inertias are given in table II). 



The body axis velocities are xelated to the carriage velocities through 
the same coordinate transformation outlined shove. All of the velocity 
components r e a lised at each step of the transformation, along with the 
expressions for the body vein velocities In terns of the space axis 
velocity, in coaplsta and snail angle linearised fom, art praaastad 


I 
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in figure 34. Again note, for the present study, 0 is used in place 
of £• Using small angle perturbations, the most significant velocity 
expression is 

v « v f • U c t (l8'< 

f 

which will be used to convert the three -degree -of -freedom equations to 
the gimbal frame of reference. 

The linearized, small perturbation lateral-directional equations 
of motion are: 

Y v v + g0 - m(v + U 0f i) - 0 

LyV + L^-jf+Ltf-O (19) 

N v v + Ify# ♦ Hfi - V - 0 

For the present study, though, these equations were modified to the 
gimbal axis notation by: first, accounting for the inertia differences 

in the rolling and yawing moment equations; second, substituting the 
expression for the body sideward velocity term. These substitutions 
yield: 

I Tf Tf - Uoflvft - «Tf + I# * 0 

*r f .f - Uo jht* ♦ 


t 
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Note: The t coefficients can be expressed as Y^, L^, and N^, 

respectively,, in the above equations. 

One additional correction must be ;.iade to account for the "non- 
lifted" mass of the model support and gimbal system. Defining as 
the total traveling mass, equal to the sum of model and linkage equip- 
ment, and m' as the total mass lifted (i.e., - the term used in 

% 

the lateral velocity equation may be modified by th-i ratio jpr. Ibis 
term effectively reduces the lateral acceleration produced by the 
given thrust, since this is less than the tot' .1 weight of the moving 
system. 

8lnoe no dynamic tests were made with control Inputs, these terms 
are not Included. Other assumptions Implicit in this form are: the 

vehicle Is In level flight; the time rate of change of the product of 
Inertia and moment of Inertia terms Is negligible. Although the effect 
of the product of inertia (lx x ) could be neglected, it was retained 
because the method used to simulate the descent angle caused the pro- 
duct of inertia to be more significant at the steeper descents. 

She space axis equations of motion in operator form, with the 
mass correction term are: 

(*Tf - jjr * d* ♦ Tqt ■ 0 

w f w # ♦ (LjB - B 2 )* ♦ 146 - e 2 )t - 0 (a) 

V* «“)* ♦ Oh ♦ W - •*>♦ • 0 
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The nontrivial solutions for Vf, 0, 'k are realized through the 
characteristic equation 



The natural modes of the aircraft motion are determined by the 
roots of the characteristic equation, that is, the values of "S" 
satisfying the above expression. Use of a space-fixed axis for the 
lateral -directional equations of motion results in a fifth-order 
system, or a fifth-degree characteristic equation. Since the con- 
stant term in the equation is equal to zero (ref. 19), one root is 
zero and remains so whenever yaw angle rather than yaw rate is used 
as a variable. 

For the restricted degree -of -freedom tests, the reduced set of 
equations that apply are: 

(l) Single degree of freedom in ♦, (v f ^ * 0) 

Ny + NjS - S 2 - 0 (23) 


^*1 tHaty ■„ < v - - * Jfe ^ ,. r 


50 


t 

(2) Ttfo degrees -of -freedom 


(a) 0 - \|r, (v f = 0) 
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DATA ANALYSIS FOR CASE 1 (LEVEL FLIGHT ) 

Single Degree -Of -Freedom 

For each test case, the above equations vere used to experimen- 
tally determine the lateral -directional stability derivatives of 
interest. Initially, for the level case, a yaw degree -of -freedom run 
was made with the model motors off (run 92, fig. 12). From the period 
of the motion = 2.75 l/sec was determined, and knowing the mech- 

anical spring constant to be - 31.7 ft-lb/rad, the yaw moment of inertia 

p 

was computed as 4.1 slug-ft . Using the average values for the period 
and damping ratio, from table III, the single-degree-of -freedom root 
was determined as S = -0.33+J2.86. Substituting into the single 
degree equation, the results are: 

= - 0.67 l/sec 

(mechanical plus aerodynamic) = - 8 . 4 l l/sec^ 

Eliminating the mechanical spring contribution, the (aerodynamic) 

is Ity * -8.41 - (- 7 . 58 ) - -O .83 l/sec 2 . 

Two Degrees -Of -Freedom 

With the single-degree-of -freedom data analyzed, the roll and 
yaw coupling was next examined using the concept of rotating time 
vectors and mode ratios as discussed In reference 27. For the parti- 
cular set of equations used, it is more convenient to use the mode 
ratios rather than solving the system characteristic equation. 
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The idea that the oscillatory motion of a system described by a 
linear differential equation can be represented by a vector rotating 
about its tail forms the basis for use of time vectors. The length 
of the vector is proportional to amplitude of the oscillatory motion 
of the system and the damped natural frequency is represented by the 
angular velocity of the vector. If the oscillatory motion Is damped, 
the vector length will decrease with time; if undamped, the vector 
length will increase with time. For multiple-degree-of -freedom 
systems, represented by a set of differential equations, the charac- 
teristic equation of the system contains an oscillatory pair of roots 
thus making it possible to represent the oscillatory mode in each 
variable by its own rotating time vector. The time vectors for the 
different variables in a particular mode will maintain a fixed-phase 
relationship with each other and rotate at the same frequency. The 
amplitude ratio and the phase angle between two variables are constant 
for a given linear system and do not depend upon the input or disturb- 
ance. The complex number relating the amplitudes of the two variables 
and the phase angle between them is called the mode ratio. Note, such 
an approach is also valid when only one mode is present in the response; 
therefore, for the analysis following, it is assumed that other modes 
of motion have damped out. 

The relationships for the mode ratio of roll angle to yaw angle, 
for the two-degree -of -freedom motion in roll and yaw, as obtained from 
equations ( 24 ) are 
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j. • 

f L t + Hs - j 


xz s 2 


S 2 - L03 


(26) 


0 S 2 - NjS - N. 

* 


N0S - ^5 S2 


(27) 


where "S" is the root of the characteristic equation corresponding 
to the mode of interest. Since the mode ratio is a complex number, 
two stability derivatives can be evaluated from each of the above 
equations. The phase difference between the roll angle response and 
the yaw angle response was -115° (i.e., 0 lags ^ by 115°), as 

determined from the data. The average amplitude ratio was 1.48, thus 

& 

resulting the mode ratio (j- = -0.6$ - jl. 34. The period and the 

damping ratio of the oscillatory mode were 4.69 seconds and 0.129, 

respectively, (S = -0.17 + jl.36). Using the values, lu = -0.4 l/sec 
P X2 

and — y * -O.O 67, the model space axis values for Ll and were 
*x * 

calculated using equation (26), and found to be: 


m. 


= +I.72 per second 
* +2.20 per second squared 

Having the values of and from the single-degree-of -freedom 
p xz 

data, and — * -0.049, either the real or Imaginary part of equa- 
ls 

tion (27) should have resulted In the same computed value for 
However, there was a slight difference, probably resulting from an 
unbalanced pitching moment because of the model mounting used as 
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discussed in reference 19 . Therefore, the stability derivative 
was added to the mode ratio, with the modification yielding: 


0 S 2 - ItyS - 


♦ 


tty + ttyS 



(27a) 


tty 

S, 


and Ity were then determined using the values of 

and — previously given. The results are: 

* 


K}, %, 



N0 = -O.58 per second 

N0 = +O.O36 per second squared 


The value of Ity is an order of magnitude less than that found in 
reference 19, thus corresponding to a very low, almost insignificant. 


untrinned pitching moment* 

The mode ratio relating the roll angle to the side velocity was 


used in the same fashion as £ described above, to compute the 

¥ 


remaining derivatives of interest. Prom the 0 - Vf degree -of -freedom 

(28) 


t 

v* 


^ S - Y V 
m* f 


ft 


A value of -4 - 2.16 was determined from the static data (Z * 27*5 

a 1 

mt « 1.83 slugs). The phase difference between the roll angle and the 
side velocity response Is 79° (*•«»> v f ***• 0 
from the data* 
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The average amplitude ratio was 0.105, resulting in the mode "atio 


0 


— = 0.02 + jO.l. The average period oi‘ the oscillatory moae was 6.1 
seconds with a damping ratio £ = -0.14, the minus sign indicating an 
unstable oscillation (S = + 0.21 + jl.08). The side force due to side 
velocity was then computed from equation (28) with the result, 

Y v ^ = -0.2 per second. Using the derivatives already determined and 
the fact that the side velocity and heading are related through the 
forward velocity, the additional stability derivatives of interest 
were computed as follows: 


Y^ =-U 0 ^Y v ^ = +3*66 ft per sec^ 



= -0.12 per ft-sec 

U °f 


N v = - = +0.045 per ft-sec 

f U 0 £ 


(29) 


Hie model space axis results for case 1, along with those for the 
additional test cases of this study, are shown in table V. 

The accuracy of the linear theory used to compute the stability 
derivatives was examined by further application of the time vector 
method described in reference 27. This w&s accomplished by constructing 
the vector polygons for various degrees -of -freedom to establish whether 
or not the polygons closed as required. Specifically, the terns of the 
particular equations used were added vectorally In accordance with the 
applicable magnitude and phase relationship. 

Figure 33 shoes this magnitude and phase relationship b etw ee n the 
roll end yew eagles end their higher order terne of Interest, the 
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resulting magnitudes for each term, shown in this figure, were then 
multiplied by the appropriate stability derivative, as dictated by 
the equations of motion, and summed vectorally. The vector polygons 
for the yaw moment equation and roll moment- equation for case 1 are 
shown in figures 56 and 57, respectively. As noted, the polygons 
close quite well indicating that the values computed and the technique 
used are credible. Similar diagrams were made for other degrees -of - 

freedom and for all of the cases of this investigation. Since all ] 

1 

polygons closed with apparently the same accuracy, only the 0 - 
relationship for case 1 has been presented as an illustration. 

Three Degrees -Of -Freedom 

No stability derivatives were computed using the available three- 
degree -of -freedom data. The purpose of presenting these data was to 
indicate the unrestrained motion of this configuration for the various 
test cases and to provide an additional means of verifying the computed 
results. The period and damping ratio for each case was determined 
from this data. For example, from figure 16 (case l), the period of the 
oscillatory motion is about 4.4 seconds, and the damping ratio is 
seen to be very low, approximately £ ■ 0.04. The roots resulting 
from these terms are S ** -0.06+41.1*. 

The 8 pace axis stability derivatives of the model, as computed 
from the data (table V), were used to solve the characteristic 
equation representing the global axes as given in equation (22). 

Since the constant tern, vu" approximately zero, for nost cases, the 
equation was simplified to a quart lc. The results for each of the cases 
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are shown in root form, being 

Case Roots (model space axes) 

1 -0.014, -1.488, +0.076+ j 1.362 

2 + 0 . 134 , -1.367, -0.u^o+jx.o;?3 

3 +0.458, -1.339, -0.090+4 1.648 

4 +0.395, -1.300, -0.042+41.663 

5 -0.028, -1.340, +0.154+41.697 

As can be seen, the value of the oscillatory root determined for case 1 , 
from the data, is in good agreement with the tabulated results. The 
three -degree -of -freedom data for the remaining cases were also examined 
and found to correlate with the solutions of the characteristic equations. 
The general characteristic equation ( 22 ) was also modified to represent 
the body axis equations of motion by simply deleting the attitude terms. 
The roots for the body axis dynamic motions of the model were then 
computed using the unadjusted body axis stability derivatives (those 
computed according to appendix II). These roots, given below, were 
used to compute the period, time to double or one-half amplitude, and 
the damping ratio associated with each mode and are presented for each 
case in tables VI to X. 

Case Roots (model body axes) 

1 - 0 . 024 , -1.280, -0.023±41.200 

2 +0.221, -1.457, -O.O27+4I.I68 

3 + 0 . 668 , - 1 . 637 , - 0 . 045 + 41.239 

4 + 0 . 587 , - 1 . 522 , - 0 . 032 + 41.265 

5 + 0 . 105 . - 1 . 531 , + 0 . 183 + 41.200 


■ * 'WW'.'v— ■ rnfmn, 


I 
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APPENDIX II 
FULL-SCALE CONVERSION 


The full-scale stability derivutxvcs calculated lor the various 
cases of this experiment were determined from the model values in the 
following manner: First, the expressions relating the space or gimbal 

axis moments to the model body axes were solved tu uiuu&i'er the model 

space derivatives to the model body axis values. Second, the model body 
axis derivatives were then converted to the full-scale values. The 
details relative to these steps will now be described. 


CONVERSION TO MODEL BODY AXES 

To determine the values of the model body axes derivatives, the 
rolling and yawing moments can be expressed as functions of the three 
aerodynamic parameters of interest during this study- That is, the 
moment dependency on roll rate, yaw rate, and sideward velocity are: 


„ , dN n . dN 

“b ^ p+ 5? 


r ^ 


( 30 ) 


The equations for p, r, and v, as shown in appendix I, can be 
substituted into the above moment expressions, resulting in: 


lg * ^ (0 cos 0 - ♦ sin 0) + (♦ cos 0 + 0 sin 0) + ~ (vf - Uo f t) 

% “ (0 cos 0 - ♦ sin 0) + ~ (* cos 0+0 oin C) + (vf - U D ♦) 

op or ov f 


( 31 ) 
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These body moment expressions can then be used in the equations relating 
the gimbal and body moments (eqs. (12) and (17) from appendix l). For 
the roll and yaw equations the relationships then become 


I t(i + P = Cl cos 2 e + L sin 9 cos e + N cos e sin 0 + N Sin "e]# 
x r xz L p r P 

+ [-Lp sin 0 cos 0 + I* co^0 - Np sir?© + N r cos 0 sin 0J* 
+ jjjv cos o + My sin + Q-Uo^Ly cos 6 " Uof^v sin 0^+ 

(32) 


2n: 


I z «i|f + P xz 0 = [K p cos 2 0 + N r sin 0 cos 0 - Lp sin 0 cos 0 - I* sin 0] 
+ Q-Np sin 0 cos 0 + N r cos 2 0 + Lp sin 2 0 - I^sin 0 cos 
+ Tn v cos 0 - Ly sin tfjvf + Q^Of^v 008 ® + ^o^^v 8 * n 


Again note: ■ -U 0 ^Ly 

»* * - u o f “v 

Since the general equations, as written from equations 21 of appendix I* 
are 


Ix'i + ^xt* • * + l4t + Urf * f ♦ 

♦ ■▼fVf + v 


( 33 ) 


the a ppropriate expressions relating the gUfeal and body derivatives ere 
realised by sinply setting the corresponding coefficients equal. Hint 1st 


% 
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2 2 
L0 = Lp cos 0 + (L r + N p )sin 0 cos 0 + N r sin 0 

2 2 

= (-Lp + N r )cos 0 sin 0 + Lj. cos 0 - Np sin 0 

N^j = Np cos^O + (N r - Lp)ain 0 cos 0 - Ly sin 2 0 

• 2 2 
Nf = -(Np + LjOsin 0 cos 0 + N r cos 0 + Lp sin 0 

L = Ly cos 0 +• N v sin 0 
v f 

N v ^ » N v cos 0 - Ly sin 0 

# -U 0 ^Lv cos 0 • Uo^v sin 0 

• -UofNv cos 0 + UqjLv sin 0 


(54) 

(55) 

( 56 ) 

( 57 ) 

( 58 ) 

( 59 ) 

( 40 ) 

( 41 ) 


The first four expressions have, as unknowns, the four body axis 
derivatives shown on the right-hand side of these equations. The 
glUbal axis values, the left-hand aide of each equation, have been 
determined from the data presented. The last two groupings of 
equations (i.e., Iy^ and I ¥ ^, or Lf and If) axe identical, thus, 
either set represents two equations with two unknowns. In addition 
to these aowent expressions, the side force equation can be ewrtnert 


in the 


rj the correspondence was slaply one to one, b ec a use 


fr 
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the side forcer due to any of the angular rates *rere considered 
negligible . 

Using the first four equations ((34) through (37) )> the body 
axis stability derivatives are solved a, 3 functions of the remaining 
body axis values. The first and fourth equations result in an expres- 
sion for Lp, while the second and third yield an expression for L r . 
These expressions, when substituted back into the general equations 
from which they were determined, result in two equations having only 
N r and Np as the unknowns. Solving the two final equations simul- 
taneously yields the body derivatives, ll r and N^, which are sub- 
sequently used to determine Lp and L r - Tor case 1, the body values 
are equal to the global values because both axes are alined when the 
pitch attitude is zero. The equations also lead to this result when 
9 • 0, thus establishing a check on the equations developed for the 
analysis. As an example, the computations for case 2 will be shown. 
Tor this esse. 


0 - 10° Iyj » -O.55 per sec 
sin 0 - 0.17 14 ■ el. 25 per sec 
cos 0 • O.965 ■ -O.55 per sec 
sin 2 0 ■ 0.03 ■ -0.55 pv sec 
cos 2 © - 0.965 

The expressions for Lp and from equations (54) to (57) ere 


• ■: )r_ « a- hk P 
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Ip * 1 $ + N^, - N r - -i.08 - N r 

(* 

h,. - L* - + Np - +1.70 + Hp 

Substituting Lp and L r Into equations (34) and (35)> the pair of 


equations resulting are 


0.34 Np - 0 > a r *» 0.10 

0.935 Hp + 0.34 - -0.65 


Solving simultaneously, yields 


N r » -0.4 per see and Up --0.62 per sec 


■which in turn from ( 42 ) yields 


Lp = -0.68 per ~ec and L r = L.16 per sec 

Either of the last two sets ol’ equations ((38) and (39), or (40) and 
(4l)) could be used to find Ly and N v . Again, for case 1, the 
equations Indicate that the glabal and body values are equal. Letting 
case 2 provide the example once more, the simultaneous set (fn * 
eqs. ( 38 ) and (39)) 


L* f - -0.1 - L y cos 10° + Ky sin 10° 

■ 40.084 ■ Mr cos K)P - Ly tin 10° 


yields 




L y = -0.11 per ft-sec 

N = +O.O 65 per ft-sec 
v 

This overall approach, that is solving the set of four equations in 
four unknovr. , and the last simultaneous set, was used to calculate 
the model body axis stability derivatives for the i-tmaining cases. 

All results are shown in tables VI to X. 

FULL-SCALE DERIVATIVES 

The second step of this procedure leading the actual aircraft's 

stability derivatives was the scaling up from the model body axis 

values, now known, to the full-scale value. Two calculations were 

performed on each derivative to acquire the final full-scale values. 

The first accounted for the difference between the scaled-down moments 

of inertia of the full-scale vehicle and the actual model; specifically, 

the rolling moment derivatives must be increased by a factor 

L.l 

while the factor must be applied to the yawing moment derivatives. 

These adjusted model results were then converted to the full-scale 
values by use of the scale factors for dynamic model similarity, as 
listed in table XI. The body axis stability derivatives for the model, 
adjusted model, and full-scale vehicle for each case are presented in 
tables VI to X. 

In order to present the dynamics of tie full-scale aircraft, the 
body axis equations of motion were adjusted to account for the fact 
that the model was flying down the test track, simulating descent, in 
such a manner that the horizon was thought of as being rotated, rather 
than having the aircraft actually descending. Only the Bidefarce 
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equation is affected by this condition, and was rewritten to include 
the appropriate components of track velocity and gravity. The modified 
sideforce equation, as noted in reference 26, becomes: 

]£y = mQv + Ur - Wp - g(sin 0 a )^ - g(cos 9 a >«J (43) 


where : 


U = Uj cos a 
W =. Uf sin a 

0 a = aerodynamic pitch angle 

between fititious horizon 
and principal X-axis 


Use of the aircraft's body axis eliminated all terms dependent on 
attitude (i.e., Ity = Y^ » Ly = N0 = 0), the cross products of inertia 

QU. 

and permitted setting — = 1. The modified characteristic equation 

m' 

then becomes: 


Y v -S 


g cos 0 a + (Uf s. i <x)S g sin a - (Uf- cos a)S 


L 3 - S £ 
P 


N r S - S c 


from which the roots for the full scale aircraft were computed. 
Although the order of the equation is quintic, the constant terms for 
each case were approximately zero, and thus eliminated. The roots 
for the resulting quartlc characteristic equation were: 
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Case Foots (full-scale body axes) 

1 -0.014, -O.76O, -O.OO5+4O.753 

2 .C.CS2, -O.736, -O.Oi3jJO.79i 

3 +O.283, -0.779, -0. 085JJ0.748 

4 +0.262, -0.715, -O.O76+0O.783 

5 +0.072, -0.744, +0.014+J0.770 j 

> * 
These roots are graphically presented on figure 38 and were used i 

to compute the period, damping ratio and time to double or one -half 

amplitude of the various modes, for each case, as noted on tables VI 

to X. It should be noted that the ad4usted model axis roots ware 

also determined using equation (44) and used to compute the periods 

damping ratios, etc., for the model: these results compared well 

with those values scaled down from the full-scale results. 
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TABLE I.- SUMMARY OF TEST CONDITIONS 


All tests conducted at i w = 40°, = 60° and model rpm 

except where rpm = 0 as noted. 


= 4000 


Fuselage 

Collective 

Descent 

Aerodynamic 

Trim 


pitch 

pitch 

angle 

pitch 

velocity 

Run 

attitude 

P.75R 

7 

attitude 

I % 

Nos. 

0 — a 



e n a a + 7 



(deg) 

(deg) 

(deg) 

(deg) j 

(ft/sec) 





I 

17.8 

91 




► 

> 

i 

0 

92** 


Degrees 

of 

freedom 



(case 1) 13*5 


10 

(case 2) 



17.9 

204 

16 

0-t-Vf 

17.8 

213 

214 

17 

♦ 

17.7 

210 

212 

18 

0-* 

17.6 

220 

224 

19 

j 

0-Vf 

18.1* 

228 

229 

20 



* Mechanical spring in place, K* —31. 7 ft-iyVad 
** rpa » 0 










TABLE I.- Concluded. 


All tests conducted at = 40°, 5^ = 60° and model rpm = 4000 

except where rpm = 0 as noted. 

" "| 

T^clage Collective Descent Aerodynamic Trim Degrees 

pitch pitch angle pitch velocity Run Figure 0 f 
attitude 0 7 attitude U f Nos. Nos. f ree dom 

0 = a ’ i 8„ = a + 7 


20 

(case 3) 


20 

(case 5) 


(deg) (deg) j (deg) & (deg) I (ft/sec) 


13-5 ! -11 



Figure 

Nos. 

Degrees 

of 

freedom 

21 


22 

0-^ 

23 

10- v f 

24 

^SL 

1 

<- 

< 

4 



29 0-t 


30 I 0-W f 


31 0-*-v f 





































REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


TABLE III.- SUMMARY OF TRANSIENT RESPONSE DATA (MEASURED) 

FOR SINGLE DEGREE OF FREEDOM 
U only) 



Run 

^-period 

(sec) 

C 


Run 

ij/ -per i od 
( sec) 



86 

2.17 

0.110 


254 

5.80 

0.175 


87 

2.17 

.11J 

•Cuce $ 

235 

4.50 

.162 


83 

2.15 

.U6 


256 

4.80 

.158 


89 

2.20 

.116 

1 

Average!'. 

5-03 

0.165 


96 

2.20 

.112 






97 

2.20 

.118 






98 

2.18 

.122 


260 

5-20 

0.20 

Case 1 








(with 

99 

2.50 

.151 


261 

5-50 

•19 

mech. 








springs ) 

102 

2.20 

. 100 

Cano 'i 

262 

5.80 

.21 


105 

2.15 

.118 


265 

3.50 

.20 


104 

2.20 

. 105 


Averages 

5.50 

0.20 


105 

2.20 

.108 






106 

2.20 

.124 


323 

4.50 

0.18 


Averages 

2.2 

0.115 


524 

5.10 

•17 


213 

5-30 

0.220 

Case 5 

325 

5.10 

.19 

Case 2 

21JR 

5.40 

.220 


526 

5.00 

.18 


Averages 

5-35 

0.220 


Averages 

4.93 

0.18 


























TABLE IV.- Concluded. 


Case 1 


Case 2 


Case 3 


Case 4 


Cam 5 


Two degree of freedom (<fi 


^-period p 

v sec) 


Averages 


Averages 


Averages 


267 

269 

270 


Averages 


294 

295 
516 
517 

327 

328 
329 


- 0.190 


-0.140 


- 0.130 


- 0.175 

- 0.190 

- 0.180 


-0.170 


-0.190 

- 0.200 

- 0.200 


r»nv-»] 


- 0.208 

-0.195 

-0.195 


v f -period 

(sec) 

C 

5-6 

-O.I 85 

6.7 

- 0.185 

5-8 

-0.197 


6.0 

-0.110 

5-8 

-0.127 

6.7 

-o.i4l 

6.0 

-0.157 


b.j 

-0.125 

6.5 

-0.140 

6.5 

-0.150 


5.4 

- 0.135 

6.6 

-0.170 

6.0 

- 0.135 


-0.180 

-0.190 

- 0.170 


[iWMfll 


- 0.190 

- 0.195 

- 0.195 


f 
























TABLE V.- MODEL SPACE AXIS Z' ABILITY DERIVATIVES 


Case 1 

Cace 2 

Cuce 3 

Cace 4 

- 0.67 

-0.55 

-0.4l 

-0.46 

- 0.83 

-1.31 

-1.66 

-1.49 

-0.58 

-0.53 

-0.13 

- 0.29 

+0.056 

-0.15 

-0.077 

-0.023 

+0.045 

+0.084 

+0.09 

+0.08 

+1.72 

+1.25 

+1.80 

+1.76 

+2.20 

+0.82 

+3.67 

+1.48 

-0.65 

-0-53 

-0.40 

-0.40 

-0.12 

-0.10 

-).o6 

-0.06 

+5.66 

+4.71 

; 

+7.75 j 

1 

+3.13 

-0.20 

-0.25 

-0.42 

-0.17 

1 

18.3 

lB.J 

18.5 

18.4 

+0.7 

-2.3 

•3.5 

-5.95 

27.3 

25.9 

27.6 

29.2 

2.16 

2.27 

2.22 

2.05 


Parameter 


l/cec 

Ity*, l/cec 2 

(» -U 0f N v ) 

Ity, l/: tc 
Ity, l/see 2 


fy, l/sec 
iy*, l/cec 2 

(** -UofLy) 

lyj**, l/cee 
L^, l/ft-sec 
ft /sec 2 

(■ ^Of^y) 

Y v > l/eec 
U 0f , ft/sfec 
Xf, lb 
Zf» Ik 


* Derivatives exist flue to use of spec# 

**Velue Hwroactaeted for each ease* 


Case 


1 


axis qrttA 
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TABLE VII.- BODY AXIS STABILITY DERIVATIVES AND PARAMETERS (CASE ?) 


Parameter 


Ly, l/ft-se'-' 

Y v , -/see 
I x , slug-ft 2 
I z , slug-ft 2 

Oscillatory mode (Dutch roll) 


Period, sec 

Time to one-half amplitude (T 3/2) 
Damping ratio (C) 

First real mode (soiral 


Time to double amplitude (Tg), 
Second real mode (rolling mode 


Time to one-half amplitude (TW2) 


Adjusted 

model 

Full-scale 

aircraft 

-0.6l 

-0.19 

-0.94 

-0.30 

+0.099 

| 

+0.003 

+2.5? 

+0-75 

-1.36 

-0.43 

-0.22 

-0.007 

-0.25 

-0.08 

1.5 

| 

150,000 

2.7 

270,000 

2.5 

7.9 

14.7 

55.0 

0.02 

0.02 

2.5 

8.4 

0.29 

0.92 


Amnmmmi 
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TABLE VIII.- BODY AXIS STABILITY DERIVATIVES AND PARAMETERS (CASE 5) 


Parameter 


I x , slug -ft 
I z , slug-ft 2 

Oscillatory mode (Dutch roll) 

Period, sec 

Time to one -half amplitude (T^/2 ) > sec 
Damping ratio (£) 

First real mode (spiral) 

Time to double amplitude (T2), sec 
Second real mode (rolling mode) 

Time to one -half amplitude (T^^)* sec 


Model 

-Adjusted 

model 

Full-scale 

aircraft 

+0.l4 

+0.21 

+0.066 

-0.50 

-0.46 

-0.15 

-0.06 

+0.10 

+0.005 

+1.65 

+3.26 

+1.05 

-0-95 

-1.90 

-0.60 

-0.087 

-0.174 

-0.006 

- OJ '2 

-0.42 

-0.15 

3.0 

1.5 

150,000 

4.1 

1 

2.7 

270,000 

5.1 

2.6 

8.4 

15.4 

2.4 

8.1 

0.04 

0.12 

0.12 

1.05 

0.76 

2.42 

0.42 

0.28 

0.89 



TABLE IX.- BODY AXIS STABILITY 


Parameter 

N r , 1/ sec 
N p , l/sec 
N v , l/ft-sec 
Lj., l/sec 
Lp, l/sec 
by, l/ft-sec 
Y y, l/ sec 

l x , slug-ft 2 

l y, slug-ft c 

Oscillatory node (Dutch Toll) 

Period, sec 

Time to one-half amplitude sec 

Damping ratio (?) 

First real mode (spiral) 

Time to double amplitude (l^eec 
Second real mode (rolling mode) 

Time to one -half amplitude (Tj j£ p Bee 


" n 


f 




TABLE X.- BODY AXIS STABILITY DERIVATIVES AND PARAMETERS (CASE 5) 


Parameter 


N r , 

l/sec 

V 

l/sec 

v 

l/ft- 

Lr> 

l/sec 

Dp> 

l/sec 

V 

l/ft-i 

v 

l/sec 


I ;: , slug-fx 
iy, SlUg-ft 2 

Oscillatory mode (Dutch roll) 

Period, sec 

Time to double amplitude (Tp), sec 
Damping ratio (£) 

First real mode (spiral) 

Time to double amplitude (Tg) , sec 
Second real mode (rolling mode) 

Time to one-half amplitude (T^), sec 


-O.ll 

-0A5 
+0.024 
+1.16 
-0.69 
-0.125 
—0 . 46 
3-0 
4.1 


M justed 1 full-scale 
model I aircraft 


-0.167 

- 0.68 

+0.056 I +0.001 
+2.52 
-1.58 
-0.25 
-0.46 

150,000 

270,000 

8.2 
49.O 
- 0.02 









TABLE XI.- SCALE FACTORS FOR DYNAMIC MODEL SIMILARITY 


Multiply full-scale property by scale factor to obtain 

model property. 



For A = 10 

Linear dimension 

A" 1 

0.1 

Area 

A" 2 

0.01 

Volume, mass, force 

A-* 

0.001 

Moment 

A" 4 

0.0001 

Moment of inertia 

A-5 

0.00001 

Linear velocity 

A-5 

0.516 

Linear acceleration 

A 0 

1 . 

Angular velocity 

A-5 

5.16 

Angular acceleration 

A 

10. 

Time 

A-5 

0.516 

Frequency 

A-5 

5.16 

Reynolds number 

A-1.5 

0.0516 

Mach number 

A-5 

0.516 


model linear dimension 


■where X ■ 

full-scale linear dimension 
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Figure 2. - Princeton Dynamic Model Track longitudinal mount 
vith one-tenth scale dynamically similar model. 
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Figure 3.- Princeton Dynamic Model Track 
lateral/directional mount vlth one- 
tenth scale dymanlcally similar model 




scale dynamic model 




HINGE POINT 



Igure 7.- Model wing airfoil section, showing slat 
flap, and aileron arrangenent. 












tt/MC 









Figure 11.- Static test data; descent condition determination 
Ramp velocity input. P. 75 R = & a = 0, and Ap = 0. 
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Figure 13- Kyn«lc teet data; directional transient respona 
One degree-of -freedom, 0=0, P. 75 R ~ * 7 > 

and IL, = 19 *^ ft/aec. 
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Figure 20.- Ftynaaic test data; lateral/directional transient 
response, three degrees -of -freedon, P-t-Tf. e - 10 , 
P.75R " 13.5°, 7 - -5°, and U, - 18.* ft/sec. 
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Figure 34. - Transformation of linear velocities from space-fixed 
to stability axis system for tilt -wing model gimbal system. 











